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Effects of the dipolar interaction on the equilibrium morphologies of a single 
supramolecular magnetic filament in bulk 
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We present a study on the equilibrium morphologies of a single supramolecular magnetic filament in a three- 
dimensional system and their dependence on the effective strength of the magnetic dipolar interactions. The 
study is based on Langevin dynamics simulations with a bead-spring chain model of freely rotating dipoles. 
Our results show the existence of three structural regimes as the value of the dipolar couping parameter is 
increased: a coil compaction regime, a coil stretching regime and a closed chain regime in which the structures 
tend progressively to an ideal ring configuration. We also found an indication of the existence of a further 
regime of closed structures that only can be observed for long enough chains. We discuss in detail the effects 
governing each regime and the structural transition between open and closed morphologies, as well as the 
absence of the multiloop configurations found under similar conditions in two-dimensional systems. 
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I. INTRODUCTION 

Supramolecular magnetic filaments are assemblies of 
nanometer sized magnetic particles linked by polymers 
or other molecules to form a chain. These systems 
are the basis of a novel nanotechnology that com- 
bines in a unique way the interesting physical proper- 
ties of magnetic nanoparticles (MNPs) with the intrinsic 
anisotropics of one-dimensional semi-flexible chains^. 

Research on MNPs has become an extremely active 
field in recent yearsSrii. In particular, the self-assembly of 
MNPs into stable structures, mainly governed by dipo- 
lar interactions, is a topic of great interest. The zero- 
field self-assembly of magnetic colloids into dipolar rings 
— with the individual magnetic moments disposed in a 
head-to-tail arrangement along the chain — was predicted 
more than four decades ago^ and since then has been 
extensively studied by means of theoretical models and 
simulations^^— . However, the first experimental obser- 
vations of self-assembled linear chains of MNPs are more 
recent^^'^^ and, in late years, have led to the finding 
of different self- assembled structures of two-dimensional 
dispersions of dipolar chains, including fiux closure struc- 
tures like rings, necklaces or other two-dimensional multi- 
loop configurationsiir— . The direct observation of closed 
ring-like structures in three-dimensional dispersions of 
MNPs is experimentally challenging, but their existence 
has been predicted theoretically by means of different 
simulation models^i^^"— . 

On the other hand, research on magnetic filaments as 
a specific system has been more scarce to date. After the 
seminal works reported in the late 1990s for the assembly 
of micrometric particle a^^'^^ , the progress on their syn- 
thesis methods has allowed the preparation of magnetic 
filaments made of particles of very different materials and 



^) Electronic mail: psanche2@icp.uni-stutt gart.de| 
Electronic mail; |joan@ifisc.uib-csic.eS| 



characteristic sizes^"— . In addition, the filaments have 
increased their degree of flexibility, a property which is 
mainly determined by the nature of the molecular links 
between the MNPs and, to some extent, by the intensity 
of their magnetic interactions. So far, most studies on 
single or highly diluted magnetic filaments in bulk have 
been restricted to the understanding of their dynamic 
properties under the action of external magnetic fields, 
paying a special attention to their application as nanofiu- 
idic propellers and actuator o^^'^^i'^^" — . The equilibrium 
configurational properties of the filaments have been gen- 
erally disregarded. Nevertheless, the chain morphology 
has been pointed as a relevant factor for some proper- 
ties of the filaments of high interest like its thermal and 
electric conductivity or its overall coercivity^. 

As a general hypothesis, it is reasonable to expect that 
the configurations observed in self- assembled chains of 
MNPs may appear also as equilibrium morphologies of 
magnetic filaments with an adequate degree of fiexibil- 
ity. Nevertheless, it is not obvious which conditions may 
lead to the different equilibrium structures or what will 
be the behavior of the filaments under conditions which 
do not lead to the self-assembly of dipolar chains. In ad- 
dition, there is still little knowledge on the equilibrium 
structures of highly diluted self- assembled dipolar chains 
and magnetic filaments in bulk. 

In a recent work^ we introduced a simple coarse- 
grained simulation model for the study of the equilib- 
rium behavior of semifiexible magnetic filaments near 
an attractive fiat surface. We reported the effect of the 
dipolar interactions on the adsorption transition and the 
existence of different equilibrium morphologies for the 
adsorbed chain depending on the temperature and the 
magnetic dipole strength. In agreement with the most re- 
cent experimental observations we found different closed 
chain morphologies on the plane, including rings and 
two-dimensional multiloop configurations. In the present 
work we take the same model as a basis to explore the 
corresponding equilibrium morphologies of magnetic fila- 
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merits in the bulk. In particular, here we intend to quali- 
tatively determine the effects of the magnetic interactions 
on the local and global configurations of a fully flexible 
chain without the geometrical constrains imposed by the 
presence of an adsorbing surface. 

The work is organized as follows: in Section HI] we re- 
view the details of the proposed model and the simulation 
method, in Section [IIII we present and discuss the results 
of our simulations and we end with the concluding re- 
marks in Section HVl 



II. SIMULATION MODEL AND METHOD 

In order to explore the equilibrium behavior of a mag- 
netic filament in bulk we have taken just the required in- 
gredients from our previous coarse-grained model^: the 
linking potential, the dipolar interaction and the steric 
repulsion between the beads. Therefore, in our model a 
magnetic filament is represented as a bead-spring chain 
of N identical particles carrying at their centers a point 
magnetic dipole, fl, which can freely rotate in any direc- 
tion. Since the model does not impose any anisotropy 
nor explicit probability rates on the reorientation of the 
dipoles, according to our previous results^ we expect a 
large head-to-tail alignment of the dipoles with respect to 
the chain backbone for large enough values of the dipo- 
lar moment. This head-to-tail arrangement is the natu- 
ral disposition of chain-like assemblies of ferromagnetic 
particles in absence of external fields. Nevertheless, this 
behavior is also found experimentally for chains of su- 
perparamagnetic particles as a consequence of a coop- 
erative effect: the increase of the dipole reversal barri- 
ers in every bead led by the external field generated by 
its neighbors^. Therefore, we consider that this model 
might be representative — at least at a qualitative level — 
for magnetic filaments made either with superparamag- 
netic of ferromagnetic particles under the conditions that 
lead to a large head-to-tail alignment of the dipoles. 

The details of our model are the following: the steric 
repulsions between the beads are modeled by means of a 
Weeks-Chandler- Anderson potential (WCA)— : 
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where r is the distance between the centers of the inter- 
acting beads, Ui,j{r) = 4e[{a/r)^^ — (a/r)^] is the stan- 
dard Lennard- Jones potential, rcut = is the shifting 
parameter selected to make the potential repulsive and 
a is the characteristic diameter of the beads. The bonds 
between adjacent beads in the chain are formed by means 
of a finite extensible non linear elastic potential (FENE) , 
defined as: 
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with Kf = 30/(T^, and rmax = 1.5(T. The use of this 
potential implies to have the linking springs attached to 



the center of the beads, therefore neglecting the eventual 
bond stretching produced by rotations of the linked par- 
ticles. The long-range magnetic interactions are repre- 
sented by the conventional point dipole-dipole potential 



fli ■ flj 3 [fli ■ nj] [flj ■ rij] 



(3) 



where = ri — fj is the displacement vector between 
the centers of the beads i, j and fli, flj are the dipole 
moments associated to each bead. In magnetic dipo- 
lar systems it is usual to define a dimensionless dipolar 
coupling parameter. A, to represent the effective inten- 
sity of the dipolar interactions. For a system of iden- 
tical dipolar spheres, this parameter is commonly de- 
fined as the ratio between the optimum magnetic en- 
ergy of two dipoles — corresponding to a close contact in 
a head-to-tail arrangement — over the energy of the ther- 
mal fiuctuations: X — nl/ (kTb^), where is the ex- 
perimental squared dipolar moment of the spheres, k is 
the Boltzmann constant, T the experimental tempera- 
ture and b is the characteristic separation distance be- 
tween the dipoles — typically, the diameter of the beads 
or, equivalently, the distance between first-nearest neigh- 
bors in chain-like systems. We can take reduced units 
for the experimental parameters by defining fj? = /ig/ce 
and T* = kT/ce, where Ce is the characteristic exper- 
imental strength of the pair interactions in our system. 
Therefore, we get the following expression for the dipolar 
coupling parameter: 



(4) 



This definition suggests that using A as the leading pa- 
rameter in a dipolar system makes the variation of or 
T* basically equivalent from the point of view of the sam- 
pling of the parameters space. Nevertheless, one must be 
aware of the limitations of this equivalence under too ex- 
treme conditions. In particular, we expect that for values 
of A beyond a characteristic threshold or, equivalently, at 
very low temperatures, the system behavior will become 
insensitive to a further increase of this parameter as an 
indication of the onset of ground state structures. Since 
our main goal is to determine the effects of the dipolar in- 
teractions, we have chosen to sample A by taking different 
values of the squared dipolar moment within a moderate 
range, G [0,20], at a constant reduced temperature, 
T* — 1. However, the dependence of A on 6 prevents in 
our model to sample this parameter in a straightforward 
way. In particular, we have a bounded but not fixed dis- 
tance between first-nearest neighbors in the chain. We 
expect such distance to depend, at least, on either the 
value of fi^ and the local degree of head-to-tail alignment 
of the dipoles. This latter dependence makes the value 
of b unknown beforehand for any simulation we could 
attempt, so we can not set A directly as a running dipo- 
lar parameter. On the other hand, it is very convenient 
to express our results in terms of A in order to facili- 
tate the comparison with other studies. Our approach 
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to this issue has been to take as the running dipolar 
parameter of our simulations and, subsequently, to calcu- 
late the corresponding values of A by applying Equation 
^ , with b being the measured average distance between 
the first-nearest neighbors in our simulated equilibrium 
configurations. Finally, for each selected value of the 
running dipolar parameter, four different chain lengths, 
iV = {10, 25, 50, 100}, have been sampled. 

Our simulations have been performed by means of 
Langevin dynamics. The translational and rotational 
Langevin equations of motion for the center of mass of a 
given bead i and its associated dipole are, respectively: 

M^^F,~rTV,+^f, (5) 
/•^-T'.-W+Cf, (6) 

where Fi and are the total applied force and torque, M 
and / the mass and the inertia tensor of the identical par- 
ticles, Tt and Tn the translational and rotational friction 
constants and finally and ^/^ are a zero-mean gaussian 
ramdom force and torque satisfying the usual fluctuation- 
dissipation relations. The reduced time t* is given by 
t* = t\/ eel (Ma^), where t is the measured time. In gen- 
eral, the values of the mass, the inertia tensor and the 
friction constants are irrelevant for the results associated 
to an equilibrium state. For simplicity, we have taken 
cr=l,e = ee = l, M=l and the identity matrix for the 
inertia tensor in order to ensure isotropic rotations. For 
the friction constants we have taken Ft = l,Tii = 3/4, 
since these values are known to produce a fast relaxation 
to equilibrium for dipolar system a^^'^^ . 

In order to enhance the statistics of the sampling, we 
have applied the replica exchange molecular dynamics 
method (REMD)-*! to our Langevin dynamics simula- 
tions by using /i^ as the replica parameter. For this pur- 
pose, we have searched for a minimal set of values of 
the squared dipolar moment, spanning the desired range 
{fil < fil < . . . < fi'^} e [0,20], that would give us 
exchange rates above 35% for all the replicas and chain 
lengths. By means of some testing runs we found a set 
of TO = 66 unequally spaced values that satisfied such 
criterium. In summary, the simulation procedure for 
each tested chain length has been the following: first, 
we have prepared each replica i by placing in a simula- 
tion cell with open boundaries a random chain with its 
corresponding magnetic dipole of squared moment fif as- 
sociated to each particle. Then, an equilibration cycle of 
Langevin dynamics with a time step St = 0.001 has been 
run for 3 • 10® steps, enough to ensure that a stationary 
regime for the total potential energy has been reached 
at each replica. After this equilibration, we have per- 
formed a cycle of 2 • 10® further steps for measurements. 
Finally, an exchange of the configurations from adjacent 
replicas has been proposed according the following energy 
criterium: 

Pa^, = min (l, e-K('^9+c^'.(A'^)-c^"(M^)-t^<.(M^)]) , (7) 



where Ui{n'j) is the total potential energy of the con- 
figuration corresponding to replica i when /u^ is used as 
the squared dipole moment of the particles. After the 
exchanges procedure has been finished, a new equilibra- 
tion cycle using Langevin dynamics has been started. 
These equilibration-measure-exchange cycles have been 
repeated until reasonably good statistics has been ob- 
tained. REMD simulations should run for enough steps 
to allow most replicas to complete some round trips 
through all the values of the exchange parameter. In 
our case, we checked that 500 equilibration-measures- 
exchange cycles were enough to allow for no less than four 
round trips. After such equilibration procedure, we still 
extended the simulations for 300 more cycles to collect 
the measures that have been finally used in our statistics. 

As a last remark on the simulation method, the limi- 
tation of our study to a single filament has made the di- 
rect summation the method of choice for the calculation 
of the long-range dipolar interactions. The simulations 
have been carried out using the coarse-grained simulation 
package ESPResSo 3.0.3^^. 



III. RESULTS AND DISCUSSION 

In order to present our results in terms of the dipo- 
lar coupling parameter. A, as defined in equation 21 the 
first property to be determined from the equilibrated fil- 
ament configurations is the average distance between the 
first- nearest neighbors, or mean bond length (6), as a 
function of the running dipolar parameter — the squared 
dipole moment, — and the chain length, N. Figure [T] 
shows the values of (6) obtained for every selected chain 
length along with the corresponding bond length of a 
dimer — i.e., a chain formed by just two dipoles — which 
is the minimal conceivable conformation in this system. 
As one can expect, for all the cases (6) decreases smoothly 
with increasing strength of the dipole moment. On the 
other hand, the chain length has a significant impact just 
for the shortest chains at high dipole moments, for which 
the observed decrease of (6) is lower. For longer chains, 
the change of this (6) is almost independent of N. The 
correct interpretation of this effect requires a detailed 
analysis of the arrangement of the dipoles in the chain as 
/i^, and consequently A, are increased. Such discussion is 
placed in Section llll CI Finally, once the dependence of A 
on (b) has been numerically established for each selected 
value of fj,^ and N, we can express all our results in terms 
of A. 

According to previous observations of self-assembled 
dipolar chains found in ferrofluids and for a semiflexi- 
ble magnetic filament near an attractive fiat surface, it 
is reasonable to expect either a large disorder in the ori- 
entations of the dipoles or a highly ordered head-to-tail 
arrangement, mainly following the chain backbone, for 
low and high values of the dipolar coupling parameter, 
respectively. A direct inspection of the equilibrium mor- 
phologies obtained in our simulations confirms this ex- 
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7V=100, A=1.7 W=100, A=5.4 JV=100, A=9.0 



FIG. 1. Mean bond length, (b), as a function of the squared 
dipole moment, fj,^ , for different filament lengths. Error bars 
are smaller than symbols. 



pectation, as well as the existence of important struc- 
tural changes with the effective intensity of the dipolar 
interactions. Figure [2] illustrates these changes for some 
selected values of A and TV. In general, for low values of A 
the chains show a shape similar to a swollen random coil, 
with a high disorder in the orientations of the dipoles. As 
A increases, the dipoles tend effectively to align with the 
chain backbone in a head-to-tail configuration, while the 
backbone stretches and becomes smoother. At even large 
values of A, the chains adopt an irregular ring-like closed 
structure, with their ends becoming permanently in close 
contact. Finally, these ringed morphologies tend to re- 
duce their irregularities as A is further increased. Shorter 
chains apparently show less backbone irregularities at 
high values of A, getting closer to the two-dimensional 
symmetry of an ideal ring. The rest of Section IIIII is 
devoted to the formal analysis of these qualitative ob- 
servations by means of different structural parameters. 



A. Characteristic equilibrium structures 

The squared radius of gyration, Rg, and end-to-end 
distance, R^, are useful parameters for the characteriza- 
tion of the global shape of chain-like structures. Figure 
3(a) shows the variation of the scaled mean value of such 
parameters as a function of A. Some selected examples 
of t heir p roba bility distributions are also shown in Fig- 
ures [3(by] and [3(c)] Two cases of the theoretical scaling 
relation oc N'^" have been assumed for the represen- 



tation of the data in Figure 3(a) the upper and middle 
panels show the values of and scaled with the 
exponent 2i/ — 6/5. This exponent corresponds to an 
ideal self-avoiding walk, which is the expected behavior 
for this system in the limit A — > 0. The lower panel of 
Figure 3(a) shows (R^) scaled with the exponent 2i/ = 2, 





JV=100, A =22.2 N=25, A =22.2 =10, A =22.2 





that corresponds to the ideal structure that we expect 



FIG. 2. Selected snapshots of equilibrium configurations ob- 
tained for different values of the dipolar coupling parameter, 
A, and chain lengths, N. The magnetic beads are represented 
as two-color spheres, with colors indicating the orientation of 
the associated dipole. 



for A ^ 1: the perfect ring. We can observe that the 
curves tend to collapse at different regions according to 
each scaling behavior: low values of A for 2i/ = 6/5 and 
large values of A for 2i^ = 2. However, it is remarkable 
that a perfect collapse has not been reached in any case. 
This result indicates that the system is far from an ideal 
behavior, specially for the limit of large values of A. 

Independently of the scaling applied to the data, the 
results shown in Figure [3] clearly indicate the existence 
of three different regimes and a transition-like struc- 
tural change within the explored range of parameters: 
for A < 2 the chain structures tend to compact slightly 
with A, an effect that increases with the chain length. 
For 2 < A < 7 this tendency is reversed, with a re- 
markable stretching of the structures until a maximum 
of the overall extension is reached. This maximum is 
absolute for the end-to-end distance in all the cases. In- 
stead, for the radius of gyration this maximum is just 
relative except for the shortest chains. At A ~ 7 the 
chain overall extension experiences a sudden drop, with 
the end-to-end distance falling to its minimum equilib- 
rium value, corresponding to the close contact separa- 
tion between two non-bonded, head-to-tail oriented par- 
ticles. In addition, the corresponding probability distri- 
butions show the existence of notably bigger fluctuations 
of both structural parameters in this region. This impor- 
tant change corresponds to the structural closure tran- 
sition that leads to the ring-like morphologies pointed 
above. Finally, for A > 7 the end-to-end distance re- 
mains in its minimum value, showing almost a delta 
function in its probability distribution, and the radius 
of gyration tends monotonically to a plateau. The value 
of this plateau is slightly lower than the corresponding 
to an ideal ring formed by N beads of effective diame- 
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FIG. 3. (a) Scaled mean squared value of the end-to-end 
distance (upper panel) and radius of gyration (middle and 
lower panels), as a function of the dipolar coupling parameter, 
A, for different chain lengths. The scaling shown in the upper 
and middle panels corresponds to the behavior expected for 
a self-avoiding walk, (-Rg) oc N^^^. The scaling in the lower 
panel corresponds to the ideal ring behavior, (-Rg) oc A''^. (b) 
Probability histograms of (R^) /N^'^ obtained for N = 100 
and some selected values of A. (c) Corresponding probability 
histograms for {R]) /7V®/^ 



ter 6, Rl 



= feV[4sin2(7r/Ar)] - [Nh/{2'K)] 
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also remarkable that the tendency to the plateau and 
its similarity to the value expected for an ideal ring are 
significantly lower for the longest chains. As an exam- 
ple, for = 10 and A « 30.39 the mean squared radius 
of gyration measured for our equilibrium structures is 
~ 1.92, while an equivalent ideal ring would have 
-Rg. ring ~ 1-95. On the other hand, for N = 100 and 
A « 30.79 we have measured (i?^) 139.25 in front of 
the ideal value R^ ^.^^^ « 186.88. Therefore, the differ- 
ence between the measured and the ideal values for i?^ 
at A ~ 30 is of the order of 10% for iV = 10 and 25% for 
A'' = 100. This effect will be better illustrated in Section 

Unci 

The strong dependence of the closed structures with 




FIG. 4. Dependence with A of the fitted exponent, 2i/, cor- 
responding to the scaling relation for the squared radius of 
gyration (-Rg) oc N'^'^ (main plot). The inset shows how the 
values of {Rg) measured for A = (data points) deviate from 
the ideal scaling, given by 21/ — 6/5 (solid line). 



the chain length discussed above explains the deviation 
with respect to the ideal scaling behavior for large val- 

The 



ues of A that has been observed in Figure 3(a) 
deviation corresponding to small values of A is instead 
a consequence of sampling too short chains. This effect 
can be evidenced by fitting the scaling exponent 2z/ to 
the available data. In Figure |4] we have represented the 
fitted exponent as a function of the dipolar coupling pa- 
rameter. The fit has been performed by using all the 
simulated chain lengths. The result is qualitatively con- 
sistent with the theoretical expectations: we can observe 
that the value of 2;/ is effectively bounded by the val- 
ues corresponding to the Flory exponent for the three- 
dimensional swollen random coil, 2iy — 6/5, and the ideal 
ring, 2^ = 2. The local peak observed at A ~ 6 corre- 
sponds to the maximum stretching of the chain prior to 
the closure transition. However, the inset of Figure H] 
evidences that the exponent fitted for A = is sHghtly 
higher than the theoretical value due to the fact that we 
have not reached the asymptotic scaling regime in our 
sampling. An accurate estimation of the scaling expo- 
nent for the Umit A -^^ would require the simulation of 
longer chains. 

Another general effect of the dipolar interactions on 
the chains is to locally stretch the backbone as A is in- 
creased, apparently in a rather independent way with 
respect to the global structure. This magnetically driven 
local bond stretching has a similar effect than the chain 
stiffness associated to bond bending potentials^. Nev- 
ertheless, it is important to keep in mind that the non- 
local structural effects of the dipolar interactions and the 
bond stiffness are completely different. In other chain- 
like systems governed by long-ranged interactions — like, 
for example, in polyelectrolyte systems — the effects of 
the intensity of such interactions on the local structure 
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FIG. 5. Probability distributions of the cosinus bond angle 
parameter corresponding to A'^ = 10 (left) and A'' = 100 
(right) obtained for different values of the dipolar coupling 
parameter. 




FIG. 6. Average bond-dipole alignment modulus, {A), as a 
function of the dipolar coupling parameter for different chain 
lengths. Its corresponding fluctuations for low values of A, 
scaled according to the sampling statistics, are shown in the 
inset. 



of the chain are usually studied by means of the persis- 
tence length^^'^°. In our case, however, the existence of 
a closure transition led by the intra-chain long-ranged 
dipolar interactions introduces a non-local dependence 
of this parameter (for a discussion on the locality of the 
persistence length see for instance refj^ and references 
therein). Therefore, we choose to analyze a simple but 
representative local parameter: the cosine bond angle 
distribution. This is defined as the probability distri- 
bution of the vector product of adjacent unitary bond 
vectors for every position i along the chain: 



Ci — bi-i,. 



(8) 



where bi^j is the unitary vector pointing in the direction 
from the center of bead i to the center of bead j. Fig- 
ure [5] shows the distributions for the cases N = {10, 100} 
corresponding to some selected values of the dipolar cou- 
pling parameter. The results show a continuous change 
with A from an almost flat distribution, with a wide do- 
main of values for Ci , to a distribution that progressively 
approaches a delta function at the point corresponding 
to an ideal ring: 



a(A > 1) ^ (5[cos(27r/iV)] 



(9) 



The almost flat distribution observed at low values of 
A indicates an insignificant correlation between adjacent 
bond vectors. For high values of A the adjacent bond 
vectors are highly correlated, as this corresponds to a 
locally stretched backbone. In all cases Ci is unable to 
take values far below -0.5 as a direct consequence of the 
steric repulsions. 

In Sections IIIIBI and IIII CI we analyze in more detail 
the structural properties found for low and high values 
of A, respectively. 



B. Properties of open structures 

We have shown that two different regimes are found for 
open structures. In order to understand this behavior, we 
have analyzed in the first place the degree of alignment of 
the dipoles with the chain backbone. This alignment can 
be easily calculated by means of the bond-dipole align- 
ment modulus. A, defined as^: 



A 



1 



N-l 



N -2 



(10) 



where fii is a unit vector parallel to the dipolar moment 
of bead i and is again a unit vector parallel to 

the displacement vector between the centers of the beads 
i — 1 and i + We expect that for any given equilibrium 
configuration found in our system this simple parameter 
will take values from 1/2 to 1. The value 1/2 will corre- 
spond to a distribution of dipole orientations completely 
uncorrelated with the chain backbone and the value 1 
to a configuration with all dipoles perfectly aligned with 
it. Figure IH shows the average and fluctuations of this 
parameter as a function of A for every chain length: a 
transition-like behavior for (A) around A 1.5 can be 
observed, with a relatively abrupt jump from the fully 
uncorrelated state into a highly correlated one and a very 
weak dependence on the chain length. For high values of 
A, the parameter is observed to remain monotonically 
approaching 1. Interestingly, the maxima in the fluctu- 
ations of the transition-like jump corresponds quite well 
with the minima found for and that separates 
the two behaviors associated to open structures. Fur- 
thermore, the value of A at which the ordering transi- 
tion seems to take place agrees quite well with the value 
known as the condition for the self-assembly of spheric 
MNPs into dipolar head-to-tail chains, A > 3^i^iS. The 
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FIG. 7. Relative average number of correlated lateral con- 
tacts, {nc)/N, for low values of A and different chain lengths. 
Inset shows the corresponding fluctuations. See the text for 
a definition of {ric}. 



self-assembly into head-to-tail conformations should be 
possible only if the thermal fluctuations are less signifi- 
cant than the magnetic energy of two perfectly aligned 
dipoles. Several previous studies on the self-assembly of 
dipolar particles have reported a negiglible probability 
of particle aggregation for values of the dipolar coupling 
parameter strictly below A < 2, neither into head-to-tail 
chain-like structures nor into three-dimensional clusters 
of correlated dipoles^"—. However, it should be noted 
that there exists at least one favoring factor in the dipole 
ordering that is absent in the self-assembly of freely inter- 
acting dipolar particles; in our system, the chain connec- 
tivity and the close contact between neighbor dipoles is 
guaranteed independently of the thermal fluctuations. In 
particular, in our model the bond length is not fixed but 
is strictly bounded to a short distance, thus the presence 
of the bonds reduces the degrees of freedom of the dipoles 
and, consequently, their configurational entropy, an effect 
that should favor the local alignment of the dipoles. 

According to the previous observations, we conclude 
that a compaction or a stretching of the open chains is 
obtained depending on the degree of alignment of the 
dipoles with the chain backbone. In particular, disor- 
dered dipoles tend to favor the chain compaction by in- 
creasing their lateral contacts and allowing the forma- 
tion of small, weakly interacting, three-dimensional clus- 
ters that reduce the overall extension of the chain. The 
monotonous decrease of the bond length with A may help 
slightly to the overall compaction in this region, but its 
impact can be neglected in front of the variations of the 
lateral contacts. These variations can be measured by 
means of a careful count that should include dipole cor- 
relations and energy criteria in order to exclude as far as 
possible the pure thermal fluctuations effects. Therefore, 
we defined the average number of lateral contacts, {uc), 
as the average number of bead pairs, that satisfy 



the following criteria: first, we exclude from the count 
the first-nearest neighbors, |j — j| 7^ 1, and the particles 
at the chain ends, i,j ^ {1,A^}; second, the particles 
should be in "close contact", thus their separation dis- 
tance should be lower than the maximum extension of 



bond. 



< 1.5cr; third, the absolute value of the 



normalized scalar product of their respective dipolar mo- 
ments should be bigger than the corresponding to a ran- 
dom distribution of dipole orientations, -/ijl > 0.5; 
finally, the pair magnetic energy of the particles should 
be strictly negative, Uj:)ip{fij; jli, fij) < 0. The combina- 
tion of these conditions, which are similar to the energy 
and entropy criteria established for the aggregation of 
magnetic particles in ferrofluids^^ , ensure that the corre- 
sponding lateral contacts are most likely to be driven by 
attractive dipolar interactions. Figure [7] shows the rel- 
ative values of this parameter measured for each chain 
length at the region of low values of A. As we expected, 
(rtc) increases with the dipolar parameter up to a max- 
imum at around A ~ 1.5, which is the dipole ordering 
point that separates the compaction and the stretch- 
ing regimes, and finally decreases for higher values of 
A. It can also be noted that longer chains tend to have 
slightly more average lateral contacts, leading to a more 
pronounced compaction regime, in agreement with the 
behavior observed in Figure [31 



C. Closure transition and properties of closed structures 

The inspection of the configurations and the behav- 
ior of the radius of gyration and the end-to-end distance 
has evidenced the existence of a closure transition in our 
model. It is known that the dipolar energy of a chain 
of > 4 dipoles arranged into a head-to-tail closed 
ring is lower than the corresponding to a head-to-tail 
stretched arrangement of the same length^. In partic- 
ular, for A^ > 4, the decrease of the dipolar energy led 
by the added head-to-tail close contact between the chain 
ends overcomes the effect of the local misalignment intro- 
duced by the ring curvature. On the other hand, a closed 
chain structure has a lower configurational entropy than 
a open one. Therefore, the closure transition separates 
the region of configurations controlled by entropy, corre- 
sponding to the open chain structures, from the energy- 
controlled region of closed structures. 

In order to analyze in more detail the closure transi- 
tion of our equilibrium structures, we choosed the total 
magnetization of the chain as the characteristic parame- 
ter. This is simply defined as the module of the vector 
sum of the unitary dipolar moments along the chain: 



M 



N 



(11) 



The behavior of this parameter is expected to be quali- 
tatively very similar to the observed for the end-to-end 
distanced^, but has an additional advantage over the lat- 





FIG. 8. Mean total magnetization as a function of the dipolar 
coupling parameter for different chain lengths and its corre- 
sponding fluctuations (inset). 



FIG. 9. Shape parameters of the chains for every sampled 
value of A and A'^. These shape parameters are defined as 
the ratios of the average principal moments of the gyration 
tensor: L2/L1 (upper panel) and Ls/L\ (lower panel). 



ter: since the domain of values sampled by M is con- 
siderably smaller than the corresponding to it is far 
easier to obtain accurate statistics for the fluctuations of 
M. Figure |5] represents the relative mean value of this 
parameter and its fluctuations. As expected, its behav- 
ior for A > 2 is analogous to what has been observed for 
R^: it grows with A up to a maximum, from which a 
sudden drop to almost zero is found around A ^ 7. Fi- 
nally, for high values of A it remains close to zero. This 
almost zero value of M indicates that the magnetic flux 
along the chain is following a nearly closed trajectory, as 
corresponds to a closed structure of head-to-tail aligned 
dipoles. Remarkably, the fluctuations of M, shown in the 
inset of Figure [51 display very clear peaks indicating the 
transition points, Aq, which move slightly towards higher 
values of the dipolar coupling parameter as the chain 
length increases: from Ao ~ 5.5 for N — 10 to Aq ^ 6.5 
for N — 100. These closure points are consistent with 
the structural phase diagram known for dispersions of 
free MNP^Si. In particular, ring-like closed structures 
are found for values of A higher than the self-assembly 
condition, A > 2.5, at low dispersion densities. Our re- 
sults indicate that there exist an important difference in 
the flux closure behavior of free MNPs and magnetic fila- 
ments. In particular, we have shown that in our model all 
the filaments remain persistently in a closed arragement 
for values of A higher than the transition point. How- 
ever, it is well known that self- assembled structures of 
free MNPs tend to be quite heterogeneous, with a coexis- 
tence of open and closed structures in a wide range of pa- 
rameters. This coexistence is not limited to experimental 
observations — in which other more complex short ranged 
interactions, hard to control and predict, might play an 
important role in the self-assembly process — but is also 
found in simulations with other minimalistic models sim- 
ilar to ours. To our best knowledge, an abrupt closure 
transition in self- assembled structures of free MNPs, lead 



by a change in A in absence of external fields, has never 
been observed and is not expected to exist due to the 
inherent polydispersity of such self-assemblies. This sug- 
gests that magnetic filaments with defined lengths might 
have a more coherent and controllable experimental equi- 
librium behavior than directly self-assembled structures 
of MNPs, thus they might facilitate the synthesis of more 
complex magnetic nanostructures. 

We have seen that in all the cases the chains remain 
in a closed configuration for A > Aq, but we need to de- 
termine how far are these closed structures from an ideal 
ring configuration. The closeness of the chain structures 
to specific ideal configurations can be characterized by 
means of different shape parameters calculated from the 
average principal moments of the gyration tensor of the 
sampled configurations, Li > L2 > L3. Here we choose 
to analyze the ratios of the average second and third mo- 
ments to the first one, L2/L1 and L3/L1 respectively. 
Figure [9] shows the behavior of these parameters with A 
for every chain length. For A = we can observe that 
our measures are in good agreement with the ratios cor- 
responding to the shape anisotropy of a swollen random 
coil, Li : L2 : L3 w 12 : 3 : 1—. For A > these shape 
parameters provide more insights on the structures of the 
different regimes. In the coil compaction region (A < 2) 
there is a very small increase of both ratios, L2/L1 and 
L3/L1, which can be observed for all the chain lengths 
and that corresponds to an slightly enhanced structural 
isotropy led by the lateral contacts. Most of the coil 
stretching region (2 < A < 5) is characterized by a de- 
crease in L3/L1 corresponding to the increased shape 
anisotropy that is expected from the local stretching of 
the chain. However, this effect becomes less evident as 
the chain length increases and L2/L1 seems to remain 
almost constant. In order to understand this behavior, 
one should take into account that in this region the fiuc- 
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tuations of the global structure grow with A, as we have 
shown in Figures 3(c) and |3(b)] as well as with N. Since 
the fluctuations make the chains to sample a wider re- 
gion of the configurational space, they tend to hinder the 
effects of the local stretching on the average values of the 
moments. The impact of the structural fluctuations and 
its dependence on the chain length become even more ev- 
ident in the region of the closure transition (5 < A < 8). 
In this case we can observe a signiflcant increase of the 
isotropy with respect to the two first principal 
evidences L2/L1, specially for the shorter chains. The 
dependence of the third dimension on the chain length 
is instead the opposite: as the chain length increases, a 
more pronounced increase of the isotropy with respect 
to the third principal axis can be observed. As shows 
L3/L1, the closed chains found persistently for A > 8 are 
relatively flat structures that tend to be closer to a two- 
dimensional ring as A increases. Nevertheless, L2/L1 ev- 
idences that these rings are quite asymmetric, having the 
shape of an irregular ellipse. It is known that semiflexible 
closed chains tend to show elliptical shapes rather than a 
symmetric, circular one^. More interestingly, the behav- 
ior of L2/L1 suggests the existence of a further structural 
regime of closed morphologies, for values of A just above 
of the closure region, that might be only observed for long 
enough chains. In particular, the curve corresponding to 
N — 100 indicates that the two-dimensional anisotropy 
of the closed structures is growing slightly with A within 
the region 8 < A < 14. However, to test this hypothesis 
further simulations with longer chains are needed. 

In our simulations, we found no trace of two- 
dimensional multiloop structures in equilibrium. 
This result differs from the experimental observa- 
tion in assemblies of MNPs in quasi two-dimensional 



systems 



17-19.21.22.24 



and the theoretical predictions of 



magnetic filaments adsorbed onto a planar surface^*. 
Smooth two-dimensional multiloop structures can be 
slightly more energetically favorable than a single ideal 
ring, provided they keep the closed head-to-tail contact 
between the chain ends and add some favorable lateral 
contacts. On the other hand, the formation of such lat- 
eral contacts between middle points in the chain induces 
a decrease in the configurational entropy with respect 
to a single ring. The balance between the changes 
in the configurational energy and entropy associated 
to the formation of multiloop structures depends on 
the system dimensionality: the presence of a steric or 
adsorbing surface imposes a two-dimensional arragement 
of the chains, leading to a significant decrease in their 
configurational degrees of freedom. Due to this entropy 
reduction mechanism, the formation of multiloop chains 
requires a smaller change of entropy in two-dimensional 
systems than in the bulk, thus it is favoured in a broad 
range of paramenters. Nevertheless, we expect that 
multiloop configurations will be observed in bulk for 
low enough temperatures, when energetic contributions 
overcome the entropic ones. 

At the beginning of Section IIIII we provided an evi- 



dence of the significant impact of the chain length on the 
mean bond length, (6), for short chains and strong dipo- 
lar interactions (see Figure [T]) . Instead, the variation of 
(6) with A is practically indistinguishable for long chains. 
On the other hand, we have evidenced that for large val- 
ues of A the dipoles show a high degree of head-to-tail 
alignment with the chain backbone which is independent 
of N (Figure [S]). Therefore, an important contribution 
to the dependence of the bond length with N for short 
chains should correspond to the behavior of the cosi- 
nus bond angle distribution and its upper bound limit 
for A ^ 1, defined in Equation As a consequence, 
shorter chains have a higher degree of misalignment even 
for the case of an ideal ring structure. In other words, 
the closed structure imposes a higher local curvature on 
the chain backbone for the shorter chains, which leads 
to an increased misalignment of the dipoles. Neverthe- 
less, the case N — 2 shown in Figure [1] evidences that 
the misalignment between first-nearest neighbors led by 
the local curvature is not enough to fully characterize the 
complex dependence of (6) on A''. Obviously, a system of 
two dipoles is not a closed structure and has no defined 
curvature. In fact, for high values of A, this system re- 
mains as a minimal rod with both dipoles in a highly 
aligned head-to-tail orientation. All these observations 
suggest consistently that there is a significant effect of 
the dipoles beyond the range of the first-nearest neigh- 
bors on the bond length and, consequently, on the values 
of A that characterize each structural region and the clo- 
sure transition. 



IV. CONCLUDING REMARKS 

We have studied the equilibrium behavior of a sin- 
gle, fully flexible supramolecular magnetic filament under 
the influence of long ranged dipolar interactions between 
their constituent MNPs, and the anisotropics and en- 
tropic effects associated to the chain-like structure. The 
study has been carried out by means of Langevin dy- 
namics simulations with a bead-spring model of linked 
but otherwise freely rotating dipoles. 

By increasing the effective strength of the dipolar in- 
teractions at a constant temperature, we observed that 
the dipoles tend to increase their alignment with the 
chain backbone in a head-to-tail arrangement, favoring 
the formation of locally stretched structures and leading 
to an increase of the orientational correlations between 
dipoles that are spatially close in the chain sequence. 
These local effects are observed together with impor- 
tant changes in the overall morphology of the chains 
that allow to clearly distinguish three different struc- 
tural regions. For low values of the dipolar coupling 
parameter the filament structure is similar to the cor- 
responding to a three-dimensional swollen random coil, 
displaying a transition-like behavior from an initial com- 
paction towards a stretching process that is associated 
to the change in the degree of the local alignment of 
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the dipoles with the chain backbone. The compaction 
region is characterized by the formation of small clus- 
ters of dipoles that establish favorable lateral contacs. 
The stretching region appears at the characteristic point 
of the transition-like alignment of the dipoles with the 
backbone, which approximately corresponds to the con- 
dition of self-assembly of free MNPs into stable clusters, 
A 2. At around A 7 the dipoles show already a high 
degree of alignment with the backbone and the chain 
experiences a closure transition, signaling the onset of 
energy-controlled closed configurations that consist in ir- 
regular elliptic ring-like structures. In general, as the 
dipolar coupling parameter is further increased, such ir- 
regular closed structures tend to approach the morphol- 
ogy of an ideal ring. However, we found an indication 
of the existence of a further intermediate structural re- 
gion, observed within a range of values of A just above 
the closure transition, in which the rings formed by long 
enough chains would tend towards a more eccentric el- 
liptic structure. 

Our results also evidence the impact of the chain length 
on its structural behavior. In general, this impact can be 
considered a complex manifestation of the strong depen- 
dence of the conformational entropy on the chain length. 
These length-dependent entropic effects may affect signif- 
icantly the values of the dipolar coupling parameter that 
characterize each structural region and, particularly, the 
closure transition. In addition we have shown that, as 
long as some flexibility is allowed for the links between 
particles, the long-range contributions beyond the range 
of the nearest neighbor ones can significantly modify the 
estimation of the average bond length and, with it, the 
effective strength of the dipolar interactions in the sys- 
tem. 

We discussed the absence of equilibrium multiloop con- 
figurations in our three-dimensional system analogous to 
the ones found theoretically and experimentally in two- 
dimensional systems. We think that the main reason for 
this absence is the importance of entropic effects in the 
three-dimensional system that prevent the formation of 
such multiloop configurations, but we still expect it to 
appear in bulk for low enough temperatures. The strong 
dependence of the entropy on the chain length and the 
evident inadequacy of the usual estimation of the effec- 
tive dipolar strength at very low temperatures suggest 
that the characterization of the multiloop regimes might 
be a rather challenging task. 

As a last remark, despite complex interplay between 
energy and entropy, all our results indicate the existence 
of well defined structural regimes, in difference to what is 
observed in the self-assembly of free MNPs. This suggests 
that magnetic filaments can be more adequate building 
blocks for the synthesis of complex magnetic nanostruc- 
tures. 
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